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1. Introduction 

In this article, we analyze the critical function of one parameter families of 
Hamiltonians near a resonance, by an approximate renormalization group method. 
We consider the following Hamiltonians with two degrees of freedom written in actions 
A = (A 1 ,A 2 ) € R 2 and angles <p = (ipi, </? 2 ) € T 2 : 

H(A,<p) = H (A)+eV(<p). (1) 

where the integrable part Hq of the Hamiltonian is 

H (A) =«• A+i(«. A) 2 . (2) 

The vector u> = (oj, — 1) is the frequency vector of the considered torus which is located 
at ft ■ A = for Ho- For e sufficiently small and for V sufficiently smooth, the KAM 
theorem ensures the persistence of this torus for a diophantine frequency ui [jjj |^. [j[ 
(see Ref. [|J for a proof in the case of Hamiltonians (@)-(§)). When we increase the 
parameter up to a critical value, the torus breaks up. The critical function lu t— > e c (to) 
is the value of the parameter e for which the invariant torus with frequency u> breaks 
up. We notice that this function vanishes at all rational values of to. For ui close to 
P/Q satisfying a diophantine condition, Carletti and Laskar Q found numerically, 
using the frequency map analysis, the following scaling law for the standard map 
critical function : 

log s c {uj) « a Pt Q log \w - P/Q | + cp,q forwsiP/Q, (3) 

with ctp t Q = For the main resonances 0/1 and 1/1, they proved analytically that 
this exponent is equal to one. 

We derive an approximate renormalization transformation 1Z : H h- > H' acting in 
some space of Hamiltonians, in order to analyze the scaling law for different one 
parameter families of Hamiltonians (|l]). This transformation gives an approximate 
value for the critical function £ c {oj). The approximate renormalization scheme we 
use is similar to the one set up by Escande and Doveil ||, f?J (see also Refs. ||, 
Complete versions of these transformations have been defined and studied numerically 
in Refs. (l(], |ll|, |l2|, [l3| [l4|, [l{|. All these transformations are composed by a sequence 
of canonical changes of coordinates : a partial elimination of the Fourier modes of the 
Hamiltonian and a rescaling procedure (shift of the Fourier modes, rescaling of time 
and of the actions) . The aim of this renormalization is to change the scale of phase 
space such that the renormalized system describes the original one on a smaller scale 
in phase space and longer time scale. 

In the space of Hamiltonians there are two main domains for the renormalization 
map : a domain where the iterations of 7Z converge to the integrable Hamiltonian Hq , 
and another domain where the iterations diverge to infinity. The two domains are 
separated by a surface S. We assume that S is a good approximation of the critical 
surface which is the set of Hamiltonians that have a torus at the threshold of the 
break-up (see Ref. H for numerical evidence). In what follows, we will identify S with 
the critical surface. 

The approximations we perform consist in mainly two approximations : we project 
the transformation on the two main Fourier modes of the perturbation and we neglect 
higher orders in s (the amplitude of the perturbation). These approximations reduce 
the renormalization transformation 1Z to a 4 dimensional map from which we can 
derive analytical results. 
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Following Ref [g , we derive analytically the equation of the critical surface S for the 
approximate renormalization transformation, and we analyze the precise shape of this 
surface near a resonance P/Q, i.e., when the frequency ui is close to P/Q- The result 
we found is that the scaling law (|^) is satisfied. The characteristic exponent ctp^Q of 
the scaling law (^) depends on the one parameter family we consider. For the following 
forced pendulum model (Escande-Doveil Hamiltonian jTJ, Q|) 

H(p, x, t) = p 2 /2 — e (cosx + cos(x — t)) , (4) 

which shares the same symmetries as the standard map, the characteristic exponent 
is the same as the one for the standard map critical function : oipq = 1/Q for 
< P < Q. However, this exponent a^g is not the same for all one-parameter family 
critical function since we will exhibit examples of Hamiltonians with different values 
(see Eq. ( p5| ) below). 

In Sec. we derive the equivalence between the critical function of a simple 
Hamiltonian model and the critical function of the following forced pendulum models : 

H(p, x, t) — ip 2 — e(cos(aa: — bt) + cos(cx — dt)), 

where (a, b, c, d) G Z 4 , a, c > and ad — be = 1. In Sec. ||, we construct the 
renormalization transformation by defining the two generators IZc and IZs of this 
transformation. We give the explicit formulas of this transformation as a function of 
the continued fraction expansion of the frequency of the torus. In Sec. [|, we derive the 
equation of the critical surface and we analyze its shape near a resonance by computing 
the characteristic exponents of the scaling laws of the critical function. In Sec [|, wc 
analyze the scaling law for the forced pendulum models in order to see the generality 
of the scaling law and its characteristic exponent. We also discuss a conjecture on the 
link between the critical function and the Brjuno function. 



2. Critical function of various types of forced pendulum models 

We consider the following families of Hamiltonian systems : 

H = -p 2 — E\ cos(ax — bt) — e 2 cos(cx — dt), (5) 

with (a,b,c,d) G Z 4 , a,c > and ad — be = 1. In particular, we notice that for 
a = c = d = 1 and 6 = 0, the Hamiltonian (|J) with e± = £2 = e is the forced 
pendulum (^). For these models, the invariant torus with frequency uj is located at 
p = uj for the unperturbed motion e\ = £2 = 0. It is equivalent to study the torus 
located at p = of the following Hamiltonian : 

H = —p 2 + up — £1 cos(ax — bt) — £2 cos(ca; — dt). 

This Hamiltonian is also mapped into a time-independent two degrees of freedom 
Hamiltonian : 

H(A, tp) = uj ■ A + ^Af ~ £1 cos(a(^i + bip 2 ) - £2 cos(c^i + dip 2 ), (6) 

where uj = {uj, — 1). This Hamiltonian belongs to the family (Q), and it is canonically 
conjugated with the following Hamiltonian : 

H(A,<p) =uj ■ A+ -(SI ■ A) 2 - /icos<y9i - / 2 cos(^ 2 , (7) 



Scaling law for the critical function of an approximate renormalization 



4 



where uj = (uj, —1) and l~i = (1, k) for some k > 0. The exact generalized canonical 
transformation which links the two models is given by three successive steps. The first 

step is a shift of the Fourier modes (A', (p') = (f^ 1 A,Tcp), where T = f ^ ^ 

and its transposed matrix T have integer coefficients and determinant 1 in order to 
preserve the periodic structure of the angles. After the shift of the Fourier modes, the 
Hamiltonian @ becomes : 

a 2 

H'(A, ip) = (d - cw)u)' ■ A + ■ A) 2 - e\ cos^i - e 2 cos ip 2 , 

where uj' = ((auj — b)/(d—cui), —1) and fi' = (1, c/a). The second step is a rescaling of 
time by a factor l/(d—cw) : we multiply globally the Hamiltonian by a factor 1 / (d—cuj) 
such that the linear term becomes uj' ■ A. The third step is a rescaling of the actions 
by a factor A = a 2 f{d— cut) : we replace the Hamiltonian H' by \H'(A/\,tp) (this 
transformation does not change the equations of motion), with A — a 2 /{d— coj) such 
that the quadratic term is equal to (fi' • A) 2 /2. Thus the image of the Hamiltonian 
H given by Eq. (|^) by this transformation is the following one : 

H"{A, tp)=u f -A+hn'- A) 2 - a (ei cos Vl + e 2 cos ^ 2 ),(8) 
2 [a — cuj) z 

which is of the same form as H with u>' = (auj—b)/(d—cuj), k = c/a, f\ = a 2 £\/{d— cut) 2 
and /2 = a 2 Eij{d — coj) 2 . 

We consider the one-parameter family of Hamiltonians (^|) with e\ = £i — e which 
is denoted F\, and the one-parameter family of Hamiltonians (0) with fx = fi = e, 
denoted F 2 . We denote £f P (oj) the critical function of the one-parameter family iq, and 
Eh (w) the one of the family F 2 . According to the equivalence between Hamiltonians (||) 
and Hamiltonians (0), the critical functions are linked by the following law : 

Sf p (u) = a~ 2 (d - cw) 2 e h (uj'), (9) 

for uj G [b/a, d/c], where uj' > is given by 

. auj — b , . 

J = -. (10) 

coj — d 

Consequently, we emphasize the role of modular transformations in the study of 
critical functions. One already knows that the classical diophantine conditions are 
invariant under modular transformations. Furthermore, the modular properties have 
been previously considered in the modular smoothing approach |l6| ] which analyzes 
the compensations between the singular behaviors at uj and 1/uj. 



3. Renormalization operators 

In this section, we derive the approximate renormalization map for a given irrational 
frequency uj > by defining two elementary operators IZc and Us- This construction 
follows the approximate renormalization developed by MacKay in Ref. However 
the resulting renormalization transformation is different since in Ref. || , the scheme is 
based on the analysis of stability of nearby periodic orbits by looking at the residues. 
We perform TZc if oj > 1 and TZs if oj < 1. The renormalization operator IZc is exact 
(it is generalized canonical transformation without any truncature) , and TZs contains 
a projection on some relevant Fourier modes of the Hamiltonian, and we neglect the 
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high orders in the perturbation parameter. 

Each of these two elementary renormalization operators acts on the following family 
of Hamiltonians : 

H(A,<p) =u> ■ A+ -(ft • A) 2 - fx cos (pi - / 2 cos</? 2 , (11) 

which belong to the family of Hamiltonians (|j). The vector u> = (to,— 1) is the 
frequency vector of the considered torus, and the vector ft = (1, k) is the direction of 
twist in the angles. We notice that if fc/1/2 = 0, the Hamiltonian (111]) is integrable. 



3.1. Definition oflZc 

The renormalization operator TZc contains a shift of the Fourier modes and a rescaling 
of time and of the actions. The shift of the Fourier modes is constructed such that it 
exchanges the mode (1,0) and (0, 1). More precisely, we require that the new angles 
tp' satisfy : 

COS ifii = COS lf2 , 
COS ip' 2 = COS ifi . 

This is performed by the following linear canonical transformation : 

(A,<p)^(A',tp') = (CA,C<p), 
where C is the following symmetric orthogonal matrix : 




The vectors oj and ft are changed into Cu> — (— 1,oj) and Cft = (k, 1). We impose 
that the images u/ and ft' of the vectors oj and ft by the renormalization TZc satisfy 
the following normalization conditions : The second component of u/ must be equal 
to -1, and the first component of ft' must be equal to 1. We rescale the time by a 
factor lo, i.e. we multiply the Hamiltonian by 1/w, and we change the sign of phase 
space coordinates (A,ip') = (— A, — tp). Then u)' is given by u)' = (1/w, — 1), and 

n' = (i,i/fc). 

The quadratic part of the Hamiltonian becomes k 2 (ft' ■ A) 2 /(2a;). In order that this 
quadratic term is equal to (ft' ■ A) 2 /2, we rescale the actions by a factor Ac = k 2 /uj, 
i.e. we change the Hamiltonian H into \cH(A/\c,(f). 
In summary, a Hamiltonian H given by Eq. ([llj) is mapped into 

H' ( A, tp) = «' • A + X - (ft! ■ A) 2 - f[ cos V x - f 2 cos <p 2 , 

where a/ = 1/a;, kl = 1/k, and 
/{ = f 2 k 2 /u, 2 , 
f 2 = hk 2 ju?. 

The renormalization operator TZc is equivalent to the 4-dimensional map given by the 
above equations : 

(fx,h,u,k)^(fi,f 2 ,l/u,l/k). 
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In order to simplify the computations, we denote f 2 = k 2 f 2 /uj 2 . The renormalization 
operators is given by : 

fx = /2, 

h = h- 

We notice that the new frequency <J = l/u> is smaller than one. 
3.2. Definition oflZ s 

The renormalization operator IZs acts on the family of Hamiltonians ([ll]) for ui < 1. 
It contains an elimination of the mode (0, 1) of the perturbation, and a rescaling 
procedure (shift of the resonances, rescaling of time and of the actions) such that the 
image of a Hamiltonian H given by Eq. (0) is of the same general form as H and 
describes the system on a smaller scale in phase space and on a longer time scale. 
We eliminate the mode (0, 1) of the perturbation, by a near-identity canonical 
transformation. We perform a Lie transformation generated by a function F(A, ip). 
The image of a Hamiltonian H is given by : 

H' = exp(F)tf = H + {F, H) + ~{F, {F, H}} + ■ ■ • , (12) 

where { , } denotes the Poisson bracket : 

r f -, = df_ dg__dg_ df 
U,9i dip' dA dip' dA' 

and the operator F acts on H like FH — {F, H}. The generating function is chosen 
linear in the actions and of the form : 

F(A,<p) = -/ a (l + fcn- A) sin (13) 

The generating function F satisfies : 

{F,H } = f 2 cosp 2 + 0{eA 2 ) 1 

where e denotes the order of f\ and f 2 , and H is the integrable part of H given by 
Eq. (||). The image of a Hamiltonian H given by Eq. ( pd| ) is : 

H' = H - fi cos ipi - f 2 cos p 2 + {F, H } 

- A^cos^i} - / 2 {F,cos^ 2 } + ~{F, {F,H }} + 0(e 3 ). 

We neglect the order 0(eA 2 ) produced by the transformation together with the orders 
0(e 3 ). These approximations are consistent since for small e, A(t) is of order e near 
the torus with frequency to. The term —f 2 {F, cos ip 2 } is independent of the actions and 
contains the modes and (0,2) that we neglect. The term — fi{F, cos <pi} generates 
the modes (1,1) and (1,-1). We neglect the last one and we keep the next relevant 
Fourier mode —\kf\f 2 cos{pi + ip 2 ). Then after the elimination of the mode (0,1), 
the new Hamiltonian is equal to : 

H' = u ■ A + (CI ■ A) 2 /2 - h cos^! - hhf 2 cos(^i + ip 2 ). (14) 

In this elimination procedure, we have neglected the Fourier modes (1, —1) and (0, 2) 
which are of the same order e 2 as the mode (1, 1) that we kept. The justification for this 
choice is that the mode (1,1) is associated with the smallest small denominator u> ■ k 
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among these three modes : \u • (1, 1)| = 1 — u < |w • (1, — 1)| = 1 + w < \u> ■ (0, 2)| = 2. 
Thus we expect that the dynamics is mainly influenced by the mode (1,1) at the 
order e 2 of the perturbation. In phase space, this means that the periodic orbit with 
frequency vector (1, 1) is closer to the torus (with frequency u>) than the periodic orbits 
with frequency vectors (1,-1) and (0,2). 

We shift the Fourier modes according to the following linear canonical transformation : 

(A,<p)~{A', v J) = (S- 1 A,S<p), 

where 

' 1 1 
1 

such that cos^i = cos<^ 2 and cos(<pi + <^2) = cos^. The vectors u and f2 are changed 
into Su: = (ui — 1, u) and SCI = (k + 1, 1). Thus the normalization conditions on these 
vectors give the image of u> and k : 

J = (1/lo) - 1, 
k' = l/(k+l). 

We rescale time in such a way that the linear term in the actions of Ho is equal to 
u)' ■ A with oj' = (u)',— 1) : we multiply the Hamiltonian H' by a factor 1/uj. The 
quadratic term of the new Hamiltonian is equal to (k + 1) 2 ($V • A) 2 / (2w). In order to 
map this quadratic part into (fi' • A) 2 /2, we rescale the actions by a factor 

A s -(fc + l) 2 /^, 

by changing the Hamiltonian H 1 into \sH'(A/\s,<p). 

After this elimination and rescaling procedures, a Hamiltonian H is mapped into 
H" = J ■ A + l(n' • A) 2 - f[ cos Vl - f 2 cos 

where 

fc(fc + l) 2 , f 

h - 2uj2 hh, 
f2 = —r- h- 

The renormalization operator TZs is equivalent to the 4-dimensional map 

(/!,/ 2 , W, k) I ► (f[,f 2 , 1/W - 1, l/(fc + 1)), 

given by the above equations. In the variables {fi, fi — k 2 f2/uj 2 ) the renormalization 
operator is given by : 

1 _ — 2fc — •' 1 -' 2 ' 

= (r^F A - 
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3. 3. Renormalization operator 1Z 

We define the renormalization operator TZ according to the value of the frequency uj 
in the following way : 

K = K C ifw>l, 
11 = lis ifw<l. 

In order to know the sequence of operators TZq an d we look at the continued 
fraction expansion of the frequency u>. We denote the entries in the continued 
fraction of uj : 

to = [a , oi, . . .] = ■ j— . 

a " + — 

The sequence of renormalization operators is (TZs'R-c) ao ^ 1 T^s(T^sT^c) ai ^ 1 T^s ■ • ■■ We 
denote lZ a = (TZsTlc) a l T^s- This renormalization operator changes the frequency 
uj = [oq , a\ , . . .] into to' = [ai , a% , . . .] , and k = [bo, b% , . . .] into k' = [a,Q , bo , bi , . . .] . The 
explicit formulas for 7Z a arc the following ones : 

f[= Pl (uj 7 k)f?f 2 , 
/a = /O2M/1, 
uj' = X/uj — a, 
k' = l/(k + a), 

where the coefficients pi and p 2 are given by 

a-i 

fc) = 2- a (l + a/fc) + I + 1)(1 - ?cj)- 2 , (15) 

;=o 

p 2 H = (l-auj)- 2 . (16) 

For a = 1, we check that pi(u>, k) and P2(w) are equal to the coefficients of IZs- The 
relations (|T5|)-(|l^) can be proven recursively, using the fact that lZ a +i = IZs'R-c'R.a- 
Moreover, the rescaling in the actions after the renormalization lZ a is equal to 
A a = (k + a) 2 /uj, and the rescaling in time is equal to 1/uj. We notice that TZ a 
has two fixed points with k = uj = (%/a 2 + 4 — a)/2 : a trivial one /1 = f% = 0, and a 
non-trivial one fx = (pip 2 )~ 1/a , h = P2{p\P2Y X 1 a ■ 

Remark : In the construction of lZ ai we have neglected several terms, some of 
them are of order 0(e 3 ). Therefore this approximation will be valid for e small, and 
in particular in the case where the frequency is close to a rational. 

4. Critical surface and scaling behavior 

4-1. Renormalization trajectories 

First, we derive the general solution for the renormalization trajectories. We compute 
the renormalization trajectories in the variables x = (x = log fx, y = log/2). We 
denote 71 (w, k) = logpx(oj,k), 72^) = logp 2 (w), and 7 = (71,72). The action of the 
renormalization operator lZ a becomes : 



x' = N a x + *y(uj,k), 



(17) 
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a 1 
1 



By iteration, we obtain the general equation for a renormalization trajectory : 

n-2 

x n = N an _ 1 ■ ■ ■ N ao x +'y{uJn-i,k n -i)+^2 ' ' ' ^oi+iTN^O.f 18 ) 



;=o 



where 



u> n = 1/cj„_i — dn-1) 
k n = l/(fc„_i + a„-i)- 

We denote p n /qn = [ao, • ■ ■ , a n -i] the nth convergent of the frequency uj (with po = 
and go — !)■ Then 



M„ = JV OB _, • • • N ao = 



Pn 



q n -l Pn-l 

The matrix N an _ 1 ■ ■ ■ N ai+1 is equal to M n Mrjj : 



Na^ ■ ■ ■ N ai + 1 =(-1) 



/+1 



9n£>i - PnQl 



PnU+l - QnPl + 1 



In-lPl - Pn-iqi Pn-m+1 ~ <2n-lPi+l 



4-2. Critical surface 

Following Ref. we determine the critical surface that separates points for which 
x n — > — oo like —q n and from those for which x n — ► +oo like q n . Thus the critical 
surface is obtained by the condition x n /q n — > when n — > oo. The first component 
of Eq. (Ill) becomes : 



n-2 

^=x+^y+y(-l)'+ 1 



pi - —qi ) ji(ui,ki) + ( — -pi+i I 72 (wj) 

9n / \9n 



7i (w„_i,fe„_i) 



We use the fact that p n /q n — > w. Then we have : 

a; + + ^ [$71 (wj, fcj) + $+172(0;;)] = °j 



(19) 



where /?/ = (— l) i+1 (pi — w<ft). We notice that the second component of Eq. ([IJ) gives 
the same equation (|19|) for the critical surface. Following Ref. one can prove that 
the sum in Eq. ( |l9[ ) converges if and only if w satisfies the Bruno condition : 

1^ <+0 °- 



2=0 



qi 
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In order to investigate the behavior of the critical surface near a rational frequency 
P/Q = [cioj ■ ■ • we assume that uj has a large entry in its continued fraction 

expansion such that : 



UJ 



(n) _ 



[ao, ■ ■ ■ , aj-i, n, aj+i, 



We consider the behavior of Eq. (19) when n goes to infinity. In the sum, we have two 
terms : one referring to 71 and the other one to 72. The dominant term for the first 
sum is obtained when I — J and is equal to 



Pj7i(uj,kj) = (3 j log 



(kj + I) ■ ■ ■ (k,j + n) l + n/kj 
(1 - ujj f ...(!_(„_ 1) WJ )2 2™ 



The remainder of the series is of order 0(\ogn/n). The continued fraction expansion 
of ojj is equal to [n, aj+i, . . .]. Then we have luj = 1/n + 0(l/n 2 ). Furthermore, kj 
is independent of n since its continued fraction expansion is [aj-i, . . . , ao, bo, b\, . . .]. 
We use the following asymptotic expansions : 



2^ l°g(^J + 1) = n\ogn — n + O(logn), 
1=1 

n-l 

- 2 ^ log(l - /wj) = 2n + 0(\ogn). 
1=1 

These expansions are obtained by using Stirling's formula. Concerning the expansion 
of (3 j, we use the formula : 

R - ° JJ 

PJ ; j 

q.i + qj-iuj 

which follows from the fact that ojj = —(qjLO — pj)/(qj—iuj — pj_x). Since uj,j = 
1/n + 0(n~ 2 ) and qj = Q, we have : 

Pj = ± + 0{n-*). (20) 
Thus the expansion of the first sum is the following one : 

Piniuu h) = 75 (log n - log 2 + 1) + O ( 

We rewrite the second sum as : 



-2^^+ilog(l - aiuji). 



1=0 

We use the fact that 1 — aicui = wjWj+i. The dominant terms of the second sum are 
the ones for I = J ~ 1 and I = J (terms containing log n). Since ojj fa 1/n, these terms 
are of order 0(logn/n). 

In summary, near a resonance P/Q, the asymptotic form of the critical surface is : 

PI, x /logn 

x + -y = --(logn-\og2 + l) + 0' 
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In order to express this behavior in terms of the distance from u) to P/Q, we have 



from Eq. (20) 



|.-P/Q|=^ + 0(n-). 
Then we have : 

P 1 

X + Q V ~ Q l0g ~ + CQ ' 
where cq = (log 2 — 1 + 21og<5)/Q. We consider the following one-parameter family 



of Hamiltonians (11) with fx = f% = £ (denoted F2 in Sec. 0). We have the following 



scaling law for the critical function of this one parameter family : 

loge c M « ^i— log |w - P/Q\ + c PtQ {k), (21) 

where c P , Q (k) = (log 2 - 1 + 21ogQ - 2Plog(JfcQ/P))/(P + Q) for P ^ 0, and 
co,i = log 2 — 1. We notice that near the main resonance 0/1 the characteristic 
exponent of the scaling law is equal to one as it has been analytically found for the 
standard map in Ref . Q , but for a resonance P/ Q the scaling law is different from the 
one of the standard map, and is characterized by a characteristic exponent 1/(P + Q). 
Moreover we notice that the exponent ap^Q of the scaling law ([?]) is independant of 
k in that case, whereas cp t Q has a term in logfc. In the integrable limit k = 0, the 
constant Cp i g(fc) and consequently s c (u>) tend to infinity. 

In summary, we have obtained a characteristic exponent ap t Q = 1/(P + Q) of the 
scaling law (|3|) for the one-parameter family F2 by an approximate renormalization 
transformation. In what follows, we analyze the scaling law (|jj|) for other types of one- 
parameter families of Hamiltonians which are linked to the family of Hamiltonians F2 
by simple canonical transformations. 

5. Forced pendulum critical functions 

The critical function of the forced pendulum (^) is linked with the one of the one- 
parameter family of Hamiltonians (jg) accor ding to Eq. (|) : 

e/p(l - w) = e/ p (w) = lu 2 eh(1/oj - 1). (22) 

Near a resonance P/Q where < P < Q for Hamiltonians (^), the corresponding 
frequency for e# which is 1/u) — 1, is near a resonance (Q — P)/P for Hamiltonians (0). 
Thus the scaling law is the following one : 

loge /p (w)w iloglw-P/QI+Cp/Q, (23) 



where 



Cp /q = i(Iog2 - 1 + 21ogQ) + 2glogg + 2 (l - log (l - P - 
We notice that Cp/q is symmetric : Cp/q — Cx-p/q- 

The scaling law (E3I) has the same characteristic exponent 1/Q as the scaling law 
found for the standard map by Carletti and Laskar ||. This exponent is consistent 
with the fact that the standard map and the forced pendulum Hamiltonian (|J) have 
the following symmetry for the critical function : e c (l — u) = e c (uj). 
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The singular behavior of the critical function for the standard map and for the 
forced pendulum model (Q) is linked to the singular behavior of the Bruno function p7[ 
B : R \ Q — > ]R + U {00} which is defined as follows : Let x denotes a real number, (x) 
is its nearest integer (i.e. (x) — [x + 1/2], where [ • ] is the integer part), and ||x|| is the 
distance between x and (x). For an irrational frequency u>, we consider the sequences 
{8k} and {(3k} defined by 

°k = l|0fc-lll> 
k-1 



Pk-l — Y\_ 



i=0 



for k > 1 and 9 = ||w||, and we define B{oj) by 



+00 

B(w) =-^2 Pk-i log 6^. 

fc=0 



with /3_i = 1. 

For the standard map, the following conjecture is supported by numerical calcula- 
tions fjj, § : 

C : The function s c (uj)/ exp(— B(u>)) is a bounded function on [0,1]. 
A strongest statement has been conjectured in Ref. |H , by requiring that this function 
is continuous. These conjectures come from the study of the radius of convergence of 
the Lindstedt series of the standard map or complex analytic maps (see for instance, 
Refs. |l^, H)), and supported by renormalization arguments |2^, 0. For the crit- 
ical function of the standard map, this problem is not solved even if recent numerical 
results obtained in Ref. || suggest that the conjecture C holds. For the approximate 
critical function we derive for the forced pendulum (||), we expect that this conjecture 
is also true since it has the same singular behavior near rationals as the one of the 
function exp(— B(u>)). 

However, the conjecture C is not satisfied in the above statement, and this comes from 
the fact that one has to take into account the right argument for the Bruno function, 
which is not always lo. In what follows, we analyze the critical function of the forced 
pendulum (^|) in order to clarify the right choice of the argument, and to reformulate 
the conjecture C in a more appropriate form. 

The critical function of the one-parameter family of Hamiltonians (|^), denoted 
e'j (w), is linked to the critical function of the Hamiltonian ^ according to Eq. (||) : 

e' fp {uj) = a~ 2 (d - cuj) 2 e h (u'), 
for lo G [b/a, d/c], where u>' = ~(auj—b)/(cuj — d). Moreover the critical function e'j p (uj) 



is linked with the critical function £j p {lo) by the following relation (see Eq. (22)) 
e' fp {u) = a- 2 ((a -c)u + d- bfe fp {io"), 

where 

a»-b (24) 



(a — c)lu + d — b 

Near a rational frequency P/Q, the scaling law (j|) is satisfied with 

a PtQ = l/((a-c)P+(d-b)Q), (25) 
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since for the one-parameter family of Hamiltonians , this exponent is equal to the 
inverse of the sum of the numerator and the denominator of the rational frequency. 
Thus the singularities of these functions e'^ p are of a different nature than the ones of 
the Bruno function B. However, e'f p (iS) has the same singularities as exp(— B(lo")) 
where lo" is given by Eq. (p4]). The conjecture C has to be replaced by the following 
one : 

C' : The function £ c (w)/ exp(— _B(o/')) is a bounded function on [b/a,d/c]. 
Practically the nature of the singularities is not changed, and the scaling law (||) is 
general with the exponent 1/Q for a typical one-parameter family of Hamiltonians : 
The two main resonances for the Hamiltonian (Q) are b/a and d/c. They correspond 
to uj" = and ui" ~ 1 which are the two main resonances of the forced pendulum 
The characteristic exponent ctp^Q close to these main resonances is equal to 1, since 
(a — c)P + (d — b)Q = 1 for P/Q = a/b or d/c. Thus the characteristic exponent 
of the scaling law (|J) close to the main resonances is equal to one for all the 
Hamiltonians (^|) even if the main resonances are not at the same place as the ones 
of the forced pendulum model (Q). Moreover, the characteristic exponent of the 
secondary resonances is equal to 1/Q where Q if the order of the resonance in the 
parameter e in the perturbation expansion. For instance, the resonance (d + b)/(a + c) 
is of order e 2 . Its characteristic exponent is equal to 1/2 since (a — c)P+ (d— b)Q = 2 
for P/Q = (d + b)/(a + c). 
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